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We review recent advances in the physics of strongly interacting charged systems functioning in water at room tem- 
perature. We concentrate on the phenomena which go beyond the framework of mean field theories, whether linear 
Debye-Hiickel or non-linear Poisson-Boltzmann. We place major emphasis on charge inversion - a counterintuitive 
phenomenon in which a strongly charged particle, called macroion, binds so many counterions that its net charge 
changes sign. We discuss the universal theory of charge inversion based on the idea of a strongly correlated liquid 
of adsorbed counterions, similar to a Wigner crystal. This theory has a vast array of applications, particularly 
in biology and chemistry; for example, the DNA double helix in the presence of positive multivalent ions (e.g., 
polycations) acquires a net positive charge and drifts as a positive particle in electric field. This simplifies DNA 
uptake by the cell as needed for gene therapy, because the cell membrane is negatively charged. We discuss also 
the analogies of charge inversion in other fields of physics. 
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I. INTRODUCTION 

Molecular biological machinery functions in water at 
around room temperature. For a physicist, this very lim- 
ited temperature range contrasts unfavorably to the rich- 
ness of low temperature physics, where one can change 
temperature and scan vastly different energy scales in an 
orderly manner. In this Colloquium, we review the re- 
cently developed understanding of highly charged molec- 
ular systems in which Coulomb interactions are so strong 
that we are effectively in the realm of low temperature 
physics. 

To be specific, imagine a problem in which one big 
ion, called a macroion, is screened by much smaller but 
still multivalent ions with a large charge Ze each (e is 
the proton charge); for brevity, we call them Z-ions. 
A variety of macroions are of importance in chemistry 
and biology, ranging from the charged surface of mica 
or charged solid particles, to charged lipid membranes, 



colloids, DNA, actin, and even to cells and viruses. Mul- 
tivalent metal ions, charged micelles, dendrimers, short 
or long polyelectrolytes including DNA - to name but a 
few - can play the role of the screening Z-ions. 

The central idea of this Colloquium is that of correla- 
tions: due to strong interactions with the macroion sur- 
face and with each other, screening Z-ions do not po- 
sition themselves randomly in three-dimensional space, 
but form a strongly correlated liquid on the surface of the 
macroion. Moreover, in terms of short range order this 
liquid is reminiscent of a Wigner crystal, as the cartoon 
in Fig. 1 depicts. Because of its central importance we 
shall use special abbreviation SCL to denote a strongly 
correlated liquid of adsorbed Z-ions. 




FIG. 1. Strongly correlated liquid (SCL) - almost a Wigner 
crystal - of Z-ions on the oppositely charged macroion surface. 
The figure is characteristic in showing the degree to which we 
are willing to ignore the microscopic details. 



Depending on the system geometry and other circum- 
stances, correlations between screening ions may appear 
in many different ways. To create some simple images in 
the reader's mind, it is useful to begin with a few exam- 
ples. One example is that in Fig. 1, which may be the 
surface of, e.g., a latex particle screened by some com- 
pact ions. With a modest leap of the imagination, we 
can envision the same picture to represent the surface of 
a DNA double helix screened by multivalent counterions, 
such as spermine with Z = 4 (Bloomfield, 1996). Here, 
we imagine DNA as a thick and long cylinder, diameter 
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2 nm and charge — e per 1.7 nm along the cylinder, or, in 
other words, with a huge negative surface charge density 
—0.9 e/nm 2 . We study correlations between point-like 
Z-ions in Sees. Ill - V. One obvious problem with the 
model Fig. 1 is that it ignores the discreteness of charges, 
both in the macroion and in the Z-ions; chemists, for in- 
stance, hardly ever accept such continuous models. In 
Sec. XI, we address the electrostatic correlations in the 
systems with discrete charges. 

As another example, consider DNA molecules which 
screen the positive surface charge of a colloid particle. 
Obviously, DNA chains play the role of Z-ions. For very 
short DNA pieces we are back again to Fig. 1, but longer 
DNA are spaghetti-like, and in this situation, correla- 
tions mean parallel arrangement, as we see in Fig. 2. 
Theoretical treatment of this problem in Sec. VI makes 
a simplifying assumption that it may be modeled as a 
system of parallel rods. 




FIG. 2. DNA adsorbed on a positively charged surface, as 
seen in an atomic force microscopy image (Mou et al, 1995). 

Electrostatic correlations are strong not only for artifi- 
cial systems involving DNA, but also for the DNA in the 
cell (Aberts et al, 1994). In particular, to organize DNA 
in chromatin (in eukariotic cells), nature uses proteins 
having large positive charges - histones. Higher levels 
of chromatin hierarchical structure can be disassembled, 
with some of the histones released, by an elevated con- 
centration of salt. The resulting most stable lower level 
structure is a bead-on-a-string necklace, as shown in the 
Fig. 3. It is called 10 nm fiber, because the beads, called 
nucleosomes, are about that big. Each nucleosome con- 
sists of a core particle, called an octamer, with a total 
charge of about +170e. As shown in the inset in Fig. 
3, the octamer is encircled 1.8 times by a DNA having 
a charge — 292e. The electrostatic interaction energy be- 
tween the DNA molecule and the histone octamer dwarfs 
the bending energy of the DNA molecule. Strikingly, a 
simple theory discussed in Sec. VII yields a similar struc- 
ture, Fig. 4, based on purely electrostatic correlations in 
a simple model. 

Correlated screening has also many useful practical ap- 
plications. The first to be mentioned is the technology 
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of dressing of DNA with polycations (Kabanov and Ka- 
banov, 1995 and 1998), positively charged star-like poly- 
mers called dendrimers (Tang et al, 1996; Kabanov ct 
al, 2000; Evans et al, 2001), or liposomes (Radler et al, 
1997) in order to produce positive complexes with DNA. 
This facilitates gene delivery through a negative cell sur- 
face membrane (Feigner, 1997). There is also the idea 
to manufacture nano-wires by attaching positive silver 
or gold colloids to DNA (Braun et al, 1998; Keren et al, 
2001). 

We view these examples as both important and con- 
vincing enough to engage in the study of electrostatic cor- 
relations between strongly charged Z-ions. Strong cor- 
relations manifest themselves in a number of ways and 
alter dramatically the whole picture of screening. In the 
familiar Debye theory, screening somewhat reduces the 
effective value of charge as seen from a finite distance in 
the outside world. With strongly correlated ions, over- 
screening becomes possible, in which case the shielded 
macroion charge is seen from outside as having the op- 
posite sign. This counterintuitive phenomenon is called 
charge inversion. 




FIG. 3. Electron microscopy image of 10 nm chromatin fiber 
(Shao, 1999). The beads are nucleosomes. The structure of 
a nucleosome is known to about 0.2 nm resolution (Luger et 
al, 1997), and the inset shows how the DNA double helix is 
bent in the nucleosome. While the overall shape of the fiber is 
perhaps due to the sample preparation procedure, the array 
of nucleosomes on the DNA is close to periodic, similar in this 
respect to the theoretical model shown in Fig. 4. 




FIG. 4. Self-assembled complex of a negative polyelectrolyte 
molecule and many positive spheres in a necklace-like struc- 
ture. On the surface of a sphere, neighboring polyelec- 
trolyte turns are correlated similar to the rods in Fig. 8 
below. On a larger scale, charged spheres repel each other 
and form one- dimensional Wigner crystal along the polyelec- 
trolyte molecule. 
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The first experimental study related to charge inver- 
sion was reported a long time ago by De Jong (1949). 
More recently charge inverted complexes of polyelec- 
trolytes were directly observed in electrophoresis exper- 
iments (see review articles by Kabanov and Kabanov 
(1995 and 1998) and references therein). It is now under- 
stood that charge inversion is a generic phenomenon and 
it is expected in all systems in which strongly charged 
ions participate in screening. It turns out to be a natu- 
ral manifestation of correlations between screening ions. 
In recent years, the phenomenon of charge inversion 
has attracted significant attention of theorists (Ennis et 
al, 1996; Wallin and Linse, 1996; Perel and Shklovskii, 
1999; Shklovskii, 1999b; Mateescu et al, 1999; Park et 
al, 1999; Joanny, 1999; Sens and Gurovich, 1999; Netz 
and Joanny, 1999a; Netz and Joanny, 1999b; Wang et 
al, 1999; Nguyen et al, 2000a-b; Messina et al, 2000; 
Nguyen and Shklovskii, 2001a-e; Chodanowski and Stoll, 
2001; Dobrynin et al, 2001; Andelman and Joanny, 2000; 
Potemkin et al, 2001; Tanaka and Grosberg, 2001a.) 

Another equally interesting manifestation of correla- 
tions is the possibility of attraction between like charged 
macroions mediated by Z-ions (Rouzina and Bloomfield, 
1996; Gronbech-Jensen et al, 1997; Levin et al, 1999; 
Moreira and Netz, 2000b). This has implications in the 
large field of self-assembly of charged biological objects, 
ranging from RNA (Woodson, 2000; Pan et al, 1999) to 
virus heads (Gelbart et al, 2000). In Sec. IX, we show 
how this attraction of like charges competes with repul- 
sion due to inverted charges, which combine to induce 
reentrant condensation of DNA or colloids out of solu- 
tions. 

Although the subject matter of our Colloquium be- 
longs to chemical and biological physics, it has remark- 
ably many deep physical analogies in other branches of 
physics, from atomic physics to quantum Hall effect; we 
discuss this in Sec. XII. 

In what follows, we re-examine the entire concept of 
screening to include correlations. Although this requires 
stepping beyond the mean field approximation, we shall 
start with some historical remarks which also serve to 
define the notations. 



II. HISTORICAL REMARKS: MEAN FIELD THEORIES 

The history of our subject is almost one century long. 
While we are unable in this brief Colloquium to discuss 
it in any depth, we shall mention three key steps. 

In the first step, Gouy (1910) and independently Chap- 
man (1913) examined the double layer at an electrode 
surface, an intensely disputed subject at the time. Fol- 
lowing Gouy and Chapman, let us consider a massive in- 
sulating macroion. To be specific, assume that its charge 
is negative, with surface density — a. Assume further 



that the only counterions in the system are those disso- 
ciated from the surface. Since a macroion is large, the 
problem of the counterion distribution near the surface 
is one-dimensional, with both counterion concentration 
N(x) and electrostatic potential 4>{x) depending on the 
distance x from the surface. Gouy and Chapman ad- 
dressed this problem by solving the following equation: 
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Here e « 80 is the dielectric constant of water, (j) s and N s 
are the potential and counterion concentration at the sur- 
face. Formula (1) follows from the Poisson equation for 
the potential <j){x) and the assumption that the ions de- 
termining the charge density are Boltzmann-distributcd 
in the same potential. Because of this self-consistency 
assumption, this equation describes a mean field approx- 
imation. It is universally called the Poisson-Boltzmann 
equation. One boundary condition d(p/dx\ x=0 — Aira/e 
follows from the fact that the field vanishes on the other 
side of the boundary (as there is no electrolyte and there 
are no charges) . The second condition is that the concen- 
tration N{x) must be normalized to a/Ze ions per unit 
area. The exact solution of the thus formulated Gouy- 
Chapman problem reads 



N(x) = 



2 7 2 



k B Te/2ne 2 Z 
(x + Xf 



(2) 



where A is called Gouy-Chapman length, which is equal 
to 



A = kBTe/2iraZe . 



(3) 



The interpretation of Eq. (3) is interesting. If a Z-ion 
were to be alone next to the charged plane, it would be 
confined by the surface field 2ira/e to such a "height" A 
that its energy change 2-KaZeX/e is about k B T - which 
leads to the correct answer Eq. (3). Other Z-ions cancel 
the field inside the macroion but double the field in the 
electrolyte at the macroion surface (see above the bound- 
ary condition at x = 0). On the other hand, every partic- 
ular ion at every moment is higher than, roughly, half the 
other ions, whereupon it finds itself in a partially screened 
field. Exact solution of the Poisson-Boltzmann equation 
indicates that these two factors cancel each other, yield- 
ing Eq. (3). 

The next step of the screening story, well known to ev- 
ery physicist, is the theory of Debye and Hiickel (1923) 
initially developed for electrolytes - the overall neutral 
mixture of mobile ions of both signs, and now widely used 
in plasma and solid state physics. Debye and Hiickel lin- 
earized the Poisson-Boltzmann equation (1) (generalized 
by introducing the sum over ion species on the right-hand 
side). Of course, linearization can be done if the potential 
is not too strong anywhere in the system, which is often 
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the case if charges involved are small enough. Debye- 
Hiickel screening leads to exponential decay of the po- 
tential around a point-like charge Q: 



<j>(r) = 9. e -r/r s 

er 

where the Debye-Hiickel radius is given by 



= ( keT 
Ts \8irN ie 2 



1/2 



(4) 



(5) 



and Ni is the concentration of monovalent salt. 

Less well known among physicists is the fact that 
Debye-Hiickel theory ignited a heated debate: Bjerrum 
(1926) commented that / exp [e 2 /r^T] r 2 dr diverges at 
r — > and, therefore, point-like charged particles are ex- 
pected to associate in neutral pairs. The discussion led 
to the realization of the important role of short range 
repelling forces. In order to prevent association of mono- 
valent ions and formation of Bjerrum pairs, the repulsion 
forces should take over at the distance which is not much 
smaller than the so-called Bjerrum length 



£ B = e 2 /ek B T , 



(6) 



which is about 0.7 nm in water at room temperature. 
(This is why hydration layer of a few water molecules 
around each ion is essential to stabilize dissociated ions). 

The last step we mention here is relatively recent, it 
has to do with non-linear screening of cylindrical charges, 
such as the DNA double helix (Onsagcr, 1967; Manning, 
1969; Oosawa, 1971). Consider a cylinder charged to the 
linear density —rj. Since the potential is logarithmic, its 
competition with entropy is quite peculiar. Indeed, re- 
leasing counterions to some distance r requires energy 
(2eZr]/e) ln(r/a) (a being the cylinder radius), while the 
corresponding entropy gain is k B T \n.{irr 2 / no 2 ) . There- 
fore, counterions are released only as long as r] < r]z, 
where 



Vz = k B Te/eZ . 



(7) 



When a cylinder is charged in excess of — r]z, some of the 
ions remain Onsager-Manning condensed on the cylin- 
der, so that its effective net charge is equal to — rjz- To 
emphasize the importance of this subject, let us men- 
tion that the DNA double helix has a bare charge den- 
sity of about — 4.2/7i, where r]i = r\z\ z =\ ( see Frank- 
Kamenetskii et al (1987) for further DNA applications). 
Onsager-Manning condensation was more accurately jus- 
tified by Zimm and Le Bret (1983). These authors ad- 
dressed non-linear Poisson-Boltzmann equation in cylin- 
drical geometry (refining earlier works - see the collection 
of papers Katzir-Katchalsky (1971)). 

Gouy-Chapman, Debye-Hiickel, and Onsager-Manning 
theories are all of mean field type, based on the Poisson- 
Boltzmann equation. This approach works well when 



screening charges are small, Z = 1. In the case of 
strongly charged macroions and Z-ions with Z ^> 1, 
however, correlations are important and mean field the- 
ory fails. Hence it is necessary to step beyond this ap- 
proximation. This is precisely the subject matter of the 
present Colloquium. Charge inversion in this context 
should be viewed as the most obvious manifestation of 
the failure of the mean field approximation. 



III. STRONGLY CORRELATED LIQUID OF 
MULTIVALENT IONS 

To begin with, let us explain why the Gouy-Chapman 
solution (2) fails at large Z. Apart from A (3), there is a 
second length scale in the problem due to the discreteness 
of charge. It is associated with the distance between 
ions in the lateral direction, along the plane. As long 
as the system as a whole is neutral, the two-dimensional 
concentration of Z-ions is n = a/Ze, and the surface 
area per ion can be characterized by a radius R such 
that ttR 2 = l/n (see Fig. 5). Thus R = (Tin)- 1 / 2 = 
(Ze/Tro-) 1 / 2 , and hence (cf. Eq. (3) 



Z 2 e 2 /eR 
knT 



(8) 



Here T is the Coulomb coupling constant, or the inverse 
dimensionless temperature measured in the units of a 
typical interaction energy between Z-ions. A system of 
monovalent ions, Z = 1, is weakly coupled, T ~ 1, this 
is why classical mean field theory applies. By contrast, 
a system in which Z-ions have large Z is strongly cou- 
pled, and we see that R becomes larger than A. For 
example, at Z = 3 and a = 1.0 e/nm 2 we get T = 6.4, 
A ~ 0.1 nm and R ~ 1.0 nm. Clearly, mean field treat- 
ment along the lines of Poisson-Boltzmann theory fails in 
this situation, since Z-ions do not affect each other when 
they are at distances smaller than R from the plane. It 
is worth emphasizing once again that it is not only the 
linearized Debye-Hiickel theory that fails in the strong 
coupling regime, but so does also the non-linear Poisson- 
Boltzmann theory. It is the mean field approximation 
that is inapplicable because of correlations between dis- 
creet charges. 

An alternative theory appropriate for the regime T » 
1 was suggested by Perel and Shklovskii (1999). The 
main idea of this theory is that at T > 1 the screening 
atmosphere is narrowly confined at the surface (see Fig. 
1), and it should be approximated as a two-dimensional 
SCL. 

A two-dimensional liquid of classical charged particles 
on a neutralizing background, the so-called one compo- 
nent plasma, is well understood (Totsuji, 1978). At zero 
temperature, it acquires the minimal energy state of a 
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Wigner crystal, shown in the Fig. 5, in which the cor- 
relation energy per ion and the chemical potential are 
given by 

s(n) ~ -LllZV/ifc = -1.96n 1/2 Z 2 e 2 /e , (9) 

mc = = L(„) = -IMZ^/eR . (10) 

an 2 

We interpret R here as the radius of a Wigner-Seitz cell 
approximated by a disc (see Fig. 5). 

At non-zero temperature, the chemical potential of one 
component plasma can be written as fi — fiid + fJ-wc + 
5 ii. Here \iid is the chemical potential of an ideal gas at 
the same concentration. Accordingly, fiwc + Sfi part is 
entirely due to correlations. Furthermore, it turns out 
that 5/i, which is the thermal correction, is negligible for 
r» 1. Although a Wigner crystal, in terms of long range 
order, melts at T w 130, the value of 6/i is controlled 
by short range order and remains negligible as long as 
r> 1. It is in this sense that a SCL of Z-ions is similar 
to a Wigner crystal. 




• • • • 

FIG. 5. A Wigner crystal of positive Z-ions on a uniform 
background of negative surface charge. A hexagonal Wigner- 
Seitz cell and its simplified version as a disk with radius R 
are shown. 

Thus, the correlation part of the chemical potential 
can be approximated by nwc (10), which is negative 
and large: —[iwc/ksT = 1.65r ^> 1. The physics of 
a large and negative fiwc can be understood as follows. 
Pretend for a moment that the insulating macroion is 
replaced by a neutral metallic particle. In this case, 
each Z-ion creates an image charge of opposite sign in- 
side the metal. The energy of attraction to the image 
is U(x) = — (Ze) 2 /4e£, where x is the distance to the 
surface. This energy is minimal when the Z-ion is placed 
next to the surface, at a distance equal to its radius a; 
therefore, the Z-ion sticks to the surface. With this idea 
in mind, consider bringing a new Z-ion to the insulating 
macroion surface already covered by an adsorbed layer of 
Z-ions (Fig. 6). This layer behaves like a metal surface 
in the sense that the new Z-ion repels adsorbed ones, 
creating a correlation hole. In other words, it creates a 
negative image. Because of the discreteness of charges, 
the adsorbed layer is a good metal at length scales above 
R only. Accordingly, attraction to the image gets satu- 
rated at x ~ R. This is why the chemical potential of an 



ion in a Wigner crystal scales as Hwc ~ — (Ze) 2 /eR. Eq. 
(10) specifies the numerical coefficient in this expression. 
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FIG. 6. The origin of attraction of a new positive Z-ion to 
the already neutralized surface. Z-ions are shown by solid 
circles. The new Z-ion creates its negative correlation hole. 

Using the concept of images, we can now understand 
the distribution of Z-ions, N(x), near the surface. To do 
this, let us extract one Z-ion from the SCL and move it 
along the x axis. As long as x -C R, its correlation hole 
does not change, and, therefore, the Z-ion is attracted 
to the surface by the uniform electric field E — 2ira/e; 
other Z-ions do not affect this attraction in any way. 
Therefore, N(x) = N s exp(— x/X) for x <C R- Here 
N s ~ n/X is the three-dimensional concentration of Z- 
ions close to the surface plane (Eq. (1)). For x ^s> R, the 
correlation hole acts as a point-like image charge, the 
corresponding interaction energy being — Z 2 e 2 /4ex. At 
x = Z 2 e 2 /Aek B T = RqT/A = XT 2 /2, the interaction with 
the image charge drops to about fcgT, i.e., negligible, and 
the Z-ion concentration becomes 

at at ( at ( l ^Z 2 e 2 \ 

(11) 

We shall further comment on the physical meaning of A*o 
after Eq. (12). Note that the correction term — Z 2 e 2 /4ex 
to Z-ion energy, which is important in the interval R <C 
x -C Ar 2 /2, is similar to the "image" correction to the 
work function of a metal (Lang, 1973). 

The dramatic difference between the exponential de- 
cay of N(x) and the Gouy-Chapman 1/(A + x) 2 -law (2) 
is due to correlation effects. Moreira and Netz (2000a) re- 
derived these results in a more formal way and confirmed 
them by Monte-Carlo simulations. Recently, Moreira and 
Netz (2001) also showed that discreteness of the surface 
charge neglected above leads to the lateral pinning of Z- 
ions. This brings Z-ions even somewhat closer to the 
surface. General direction of this effect can be under- 
stood from the limit (although unrealistic for strongly 
charged surface) when the distance of closest approach 
between a discrete surface charge and a Z-ion is so small 
that they form isolated Bjerrum pairs (see (6)). 

At larger x, correlations and interactions with im- 
age charges are unimportant, and the Poisson-Boltzmann 
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equation applies. In this region N(x) varies so smoothly 
that N(x) — No provides an effective boundary condition 
for the Poisson-Boltzmann equation (Shklovskii, 1999b). 
At this stage we must remember that in a real physical 
situation there is always some concentration of Z-ions, 
N, in the surrounding solution. It can be either larger 
or smaller than No. In the latter case, the surface is 
overcharged, as we show in Sec. IV below. 

IV. CORRELATION-INDUCED CHARGE INVERSION 

Let us return once again to the physical argument il- 
lustrated by Fig. 6. It explains why an extra Z-ion 
may be attracted to the macroion surface despite the fact 
that it is already neutralized by the previously adsorbed 
Z-ions. What happens if another Z-ion approaches? 
Clearly, correlation effect will keep providing an attrac- 
tive force for this and subsequent Z-ions, but it will have 
to compete with the repulsive force which is simply due 
to the fact that macroion has already too many Z-ions 
adsorbed and the whole complex is, therefore, positively 
charged. Thus, the question is this: what is the equilib- 
rium amount of (over)charge? 

One useful way to think about it is to realize that cor- 
relation mechanism provides voltage to drive overcharg- 
ing, but the actual amount of (over)charge depends on 
both voltage and capacitance. Since the latter depends 
strongly on the geometry, we will have to explore several 
cases - spherical macroion, cylindrical, etc. 

Another, equivalent, view, involves the comparison of 
chemical potentials of adsorbed Z-ions and Z-ions in the 
bulk solution. This approach immediately suggests that 
in equilibrium the total charge, Q* , depends on the con- 
centration, N, of Z-ions in the surrounding bulk solution. 
Here Q* is the net charge of the entire complex which in- 
cludes bare charge of the macroion, — Q < 0, and the 
proper (determined by the equilibrium condition) num- 
ber of adsorbed Z-ions with charge Ze > each. 

Let us see now how we can implement the condition 
of equal chemical potentials for the case of a spherical 
macroion, with radius r. As regards adsorbed Z-ions, we 
argue 1 that their chemical potential is only different from 



x The argument goes as follows. Spherical surface of the 
macroion has bare surface charge density —a = —Q/4irr 2 . 
Let us pretend to place there, along with real charge —a, 
also two imaginary spheres, with uniform charge densities 
a* = Q*/47rr 2 and — a*. Their total charge is zero, so they 
have no effect. However, we can now think that the Z-ions are 
adsorbed on the sphere with the charge density — a — a* , and 
with this sphere they form a neutral SCL, quite like that con- 
sidered in Sec. HI. The remaining sphere has the charge den- 
sity a* and creates spherically symmetric field with potential 



the one-component plasma expression (10) by the energy 
of the Z-ion in the potential ru(0) = Q* jer created by the 
net charge Q* . Therefore, equilibrium condition reads 
l-Md + Mwc + Zetp(Q) = Ub, where Ub is the bulk chemical 
potential. To determine Q* from here, we first note that 
Mid — Mb — kBThi(N a /N); we further express Mwc m 
terms of A (see Eq. (11)), and finally obtain 

Q* = ^-k B T\n{N/N ) . (12) 
Ze 

Clearly, the net charge Q* is indeed positive when A > 
N , i.e. it has the sign opposite to the bare charge Q. 

The result (12) sheds also light on the meaning of the 
quantity Ao defined above, in Eq. (11): this is the 
concentration of Z-ions in the surrounding bulk solu- 
tion at which macroion is exactly neutralized by the ad- 
sorbed Z-ions. This concentration is very small because 
\^wc\/k B T > 1. For example, N = 0.3 mM and 0.8 /jM 
for Z = 3 and 4 respectively (1 M w 6 x 10 20 cnT 3 ). 
Therefore, it is easy to increase charge inversion by in- 
creasing A. How far does it go? At large enough A, 
translational entropy terms fib — Hid become negligible 
compared to Hwc, yielding 

Q*/er = i/>(0) = \ii WC \/Ze. (13) 

Expressing R and |mwc| through Q and Z with the help 
of Eq. (10) (and remembering that a* = —a + Zen), 
Shklovskii (1999b) arrived at the prediction for the max- 
imal inverted charge for the spherical macroion which can 
be achieved by increasing concentration of Z : 1 salt: 

Q* = 0.83^/QZe. (14) 

This charge is much larger than Ze, but still is smaller 
than Q because of limitations imposed by the large charg- 
ing energy. For example, for Q = lOOe, Z = 4, we get 
Q* = 17e. Eq. (14) was recently confirmed by numerical 
simulations (Messina et al, 2000; Tanaka and Grosberg, 
2001a). Further increase of charge inversion beyond the 
level dictated by Eq. (14) is achievable with the help of 
a monovalent salt (see Sec V). 

Applied in cylindrical geometry, similar arguments 
lead to a revision of the conventional Onsager-Manning 
condensation theory (Sec. II) when dealing with multi- 
valent Z-ions. Consider a cylinder with a negative lin- 
ear charge density — r\ and assume that n > rjz ■ Mean 
field Onsager-Manning theory (7) predicts rf = —r)z- 
By contrast, Perel and Shklovskii (1999) showed that a 



on the surface ij)(0). We emphasize that the macroscopic net 
charge <r* does not interact with the one component plasma, 
because the potential ^(0) is constant along the surface, while 
one component plasma is neutral. 
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correlation-induced negative chemical potential nwc re- 
sults in 



V = ~Vz 



In (4/irZ 6 N£%) 



(15) 



where £b is the Bjcrrum length (6). This result repro- 
duces the Onsager-Manning one (7) only at extremely 
small values of N, which are unrealistic at Z > 3. On 
the other hand, at N = No the net charge flips sign, re- 
sulting in charge inversion at N > iVo (which is absent 
in Onsager-Manning theory). At large enough TV, the 
inverted charge density i]* can reach fcgTe/e = rji. 



V. ENHANCEMENT OF CHARGE INVERSION BY 
MONOVALENT SALT 

Most water solutions, particularly biological ones, con- 
tain significant amounts of monovalent salt, such as 
NaCl. Correlations between these monovalent ions are 
negligible, and, therefore, their only role is to provide 
Debye-Hiickel screening with a decay length r s (5). This 
screening makes charge inversion substantially stronger. 
Indeed, screening by a monovalent salt diminishes the 
charging energy of the macroion much stronger than the 
correlation energy of Z-ions. Furthermore, in a sufficient 
concentration of salt, the macroion is screened at the 
distance smaller than its size. Then, the macroion can 
be thought of as an over-screened surface, with inverted 
charge Q* proportional to the surface area. In this sense, 
overall shape of the macroion is irrelevant, at least to a 
first approximation. Therefore, we consider here a sim- 
pler case: screening of a planar macroion surface with a 
negative surface charge density ~<r by a solution with a 
concentration N of Z : 1 salt and a large concentration 
Ni of a monovalent salt. 

Nguyen et al (2000a) calculated analytically the de- 
pendence of the charge inversion ratio, cr*/cr, on r s , 
in two limiting cases r s 3> Ro an d r s <C Ro, where 
i?o = (na/Ze)^ 1 / 2 is the radius of a Wigner-Seitz cell 
at the neutral point n — a/Ze. At r s 3> Ro calcula- 
tion starts from Eq. (13). The electrostatic potential of 
a plane with the charge density a* screened at the dis- 
tance r s reads ip(0) = 47rer*r s . At r s 3> Ro screening by 
monovalent ions does not change Eq. (10) substantially 
so that we still can use it in Eq. (13) which now describes 
charging of a plane capacitor by voltage \fiwc\/Ze. This 
gives 



a* /a = 0Al(R o /r s ) « 1 (r a » Ro) . 



(16) 



Thus, at r s » Ro inverted charge density grows with 
decreasing r s . 

Now we switch to the case of strong screening by mono- 
valent salt. To begin with, let us assume that screening is 
already so strong that r s « J)o, but energy of SCL is still 



much greater than ksT per Z-ion. In this regime, free 
energy consists of Debye-Hiickel screened nearest neigh- 
bor repulsion energies of Z-ions and the attraction energy 
of Z-ions to the charged surface: 

F = (3nZ 2 e 2 /eA) exp(-A/r s ) - 4ir<rr s Zen/e , (17) 

where A = (2/v / 3) 1 ^ 2 n^" 1 / 2 is the lattice constant of the 
hexagonal Wigner crystal (Fig. 5). Minimizing F with 
respect to n one arrives at 



a* _ it ( R a 

° ~ 2V3 \r s \n(R /r s ) 



(r s « Ro) . (18) 



Thus a* /a grows with decreasing r s and can become 
larger than 100%. At r s ~ R , Eq. (16) and Eq. (18) 
match each other. As we see, a* continues to grow with 
decreasing r s . This is because the repulsion between Z- 
ions becomes weaker, so that it is easier to pack more 
Z-ions on the surface. Of course, when r s decreases even 
further, binding energy of Z-ions becomes small, SCL 
dissolves, and charge inversion disappears. 

The results above are in a good agreement with sim- 
ulations. Terao and Nakayama (2001) reported the re- 
sults of a Monte Carlo simulation for the system consist- 
ing of a macroion with charge Q = 20e surrounded by 
20 monovalent counterions and 1500 ions of 2 : 1 elec- 
trolyte. Tanaka and Grosberg (2001a) performed molec- 
ular dynamics simulations with a spherical macroion of 
the charge Q = — 28e, spherical Z-ions (2 < Z < 7), and 
up to 500 monovalent ions; the system is neutral overall, 
which determines the number of Z-ions to be between 
180 and 52. Simulations confirmed the strong adsorp- 
tion of the overcharging amount of Z-ions on the sur- 
face of macroion. On a more quantitative level, Tanaka 
and Grosberg (2001a) examined the dependence of the 
inverted charge on the ionic strength, /, and found the 
crossover between Q* cx \fl cx l/r s and Q* cx I cx 1/r^, 
consistent with Eqs. (16) and (18). Tanaka and Gros- 
berg (2001a) attempted also to maximize the charge in- 
version ratio Q*/Q. In agreement with the theoretical 
views presented above, the growth of charge inversion is 
capped when correlations between Z-ions are suppressed. 
We have mentioned eventual reduction of correlations for 
large monovalent salt concentration, or small r s , when the 
SCL evaporates. If one tries to increase Z instead of low- 
ering r s , then, at large Z, correlations become suppressed 
because monovalent ions condense on Z-ions, forming Z'- 
ions with smaller net charge Z'. This effect is clearly 
seen in Fig. 7 (see Nguyen et al (2000b) for the conditions 
under which this phenomenon is and is not important). 
Nevertheless, Q*/Q up to about 150% is easily observed. 
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FIG. 7. The snapshot of the system simulated by Tanaka 
and Grosberg (2001a). The spherical macroion is grey, Z- 
ions are white, and monovalent negative ions are black. Bare 
charge of the spherical macroion is Q = — 28e, Z = 7, and 
r ~ 30. There are 52 Z-ions and 336 monovalent ions in the 
simulation domain. What is clearly seen is the formation of 
Z'-ions: because Z is so large, Z-ions adsorb monovalent ions. 
This reduces correlations between Z-ions and restricts charge 
inversion. Nevertheless, bare charge of the complex shown in 
the figure is +16e. 



boring rods results in the negative net surface charge den- 
sity — er* similar to Eq. (18) (Netz and Joanny, 1999a; 
Nguyen et al, 2000a): 



\n(rj/ar 8 ) 



r s <$:A . 



(19) 



Here the applicability condition involves Ao = rj/a, 
which is the distance between rods when they neutral- 
ize the plane; only at r s <C Aq is the overcharged plane 
linearly screened by monovalent salt. 

Speaking about DNA, we have already discussed in 
Sec. II that the DNA charge density, — rj, is such that 
about three quarters of it is compensated by positive 
Onsager-Manning-condensed monovalent ions. In other 
words, the net charge of DNA in the bulk solution is 
77* = —rji (7). It turns out that at r s <C A the result 
(19) applies to DNA with the only correction of replacing 
-r\ with rf = -r)i = -k B Te/e (Nguyen et al, 2000a). 



VI. SCREENING OF A CHARGED PLANE BY 
POLYELECTROLYTES 

A practically important class of Z-ions are charged 
polymers, i.e., polyelectrolytes. Let us start with a rigid 
polyelectrolyte and discuss charge inversion caused by ad- 
sorption of long rod-like Z-ions. For example, the moder- 
ately long (up to about 50 nm, or about 150 base pairs) 
DNA double helix can be well approximated as a rod. 
Actin is another example of an even more rigid poly- 
electrolyte. Apart from the uninteresting regime of ex- 
tremely small macroion surface charge density (in which 
case the elongated shape of molecules is irrelevant, ren- 
dering our previous results applicable), charged rods ad- 
sorbed at the surface tend to be parallel to each other due 
to the strong lateral repulsion. In other words, there is 
the short range order of a one-dimensional Wigner crys- 
tal with lattice constant A in the direction perpendicular 
to the rods (Fig. 8). 

To make the signs consistent with the case of DNA, we 
assume that the polyelectrolyte charge is negative and 
equal to —77 per unit length, while the macroion surface 
is a plane with positive charge density a. We assume also 
that there is a certain concentration of monovalent salt, 
Ai, in the solution, corresponding to the Debye screen- 
ing radius (5). To begin with, let us assume that charge 
density of the rods, —77, is below the Onsager-Manning 
threshold, Eq. (7), and let us apply the Debye-Huckcl 
approximation to describe the screening of the charged 
surface by the monovalent salt. We can then directly 
minimize the free energy of the one-dimensional crystal of 
negative rods on the positive surface written similarly to 
Eq. (17). Then the competition between the attraction 
of the rods to the surface and the repulsion of the neigh- 




FIG. 8. Rod-like negative Z-ions, adsorbed on a positive uni- 
formly charged plane. 

Thus the inversion ratio grows with decreasing r s as in 
the case of spherical Z-ions. At small enough r s and a, 
the inversion ratio can reach 200% before DNA molecules 
are released from the surface. It is larger than for spher- 
ical ions, because in this case, due to the large length of 
the DNA helix, the correlation energy remains large and 
the Wigner crystal-like short range order is preserved at 
smaller values of <rr s . Nguyen et al (2000a) called this 
phenomenon "giant charge inversion." 

Let us switch now to the opposite extreme of weak 
screening by a monovalent salt, r s 3> Aq. In this case, 
screening of the overcharged plane by monovalent salt 
becomes strongly nonlinear, with the Gouy-Chapman 
screening length A* = eksT / (2nea*) much smaller than 
r s . Furthermore, some of the positive monovalent ions 
Onsager-Manning condensed on DNA are released from 
it upon adsorption, as the plane repels them (Park et al, 
1999; Gelbart et al, 2000). As a result the absolute value 
of the net linear charge density of each adsorbed DNA, 
77*, becomes larger than 771. To determine a* and 77*, 
Nguyen et al (2000a) considered two equilibrium condi- 
tions, dealing with chemical potentials of rods and small 
ions, respectively. As a result, they arrived at the follow- 
ing beautiful formulae valid at r s 3> A a : 
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^ = -^cxp(-,/ln^ln-^ , (20) 

» / Info /a) A . . 

V*=Vi\ L ( A L x , r a »4>. 21 
y m(A /27ra) 

At r s ~ A /27r we get 77* ~ rji, A* ~ r s and a*/cr ~ 
r]i / (2irr s cr) so that Eq. (20) crosses over smoothly to the 
strong screening result of Eq. (19). 

So far in this Section we assumed a rod-like polyelec- 
trolyte. Let us now discuss how chain flexibility affects 
charge inversion. We argue that the results for rod-like 
Z-ions are remarkably robust. To begin with, consider 
a polyelectrolyte having several charged groups per each 
persistence length. We argue that our results remain 
valid as long as adsorption energy per one persistence 
length is large compared to fceT. Indeed, under this con- 
dition even flexible polyelectrolyte chains lay flat on the 
surface, in which case they are ordered in Wigner crystal- 
like SCL and, therefore, behave similarly to rods. 

Dobrynin et al (2001) addressed the opposite extreme, 
namely, weakly charged polyelectrolytes, with so small 
a fraction of charged monomers, /, that a link between 
two neighboring charges is already a flexible polymer; 
in other words, the distance between charges is larger 
than the persistence length. It was discovered by de 
Gennes et al (1976) that a weakly charged polyelectrolyte 
chain in a bulk solution consists of electrostatic blobs. 
Inside each blob polymer is only marginally perturbed 
by Coulomb interactions, while chain of blobs is fully 
stretched, rod-like. Dobrynin et al (2001) have found 
that this blob structure remains valid for the adsorbed 
chains, which form SCL of effective rods of blobs. Speak- 
ing of charge inversion, this means that Eq. (19) remains 
valid for the weakly charged chains, provided r\ is re- 
placed with a linear charge density of the string of blobs 
V* = (fcekBT/l 2 ) 1 / 3 , where I is the chain persistence 
length, and r s is larger than a blob size. 

With increasing a, adsorbed rods, either real or blob 
ones, start to touch each other leading to multilayer ad- 
sorption. It is only in this regime that the real and blob 
rods behave differently, as we discuss in Sec. VIII. 

VII. POLYELECTROLYTES WRAPPING AROUND 
CHARGED PARTICLES 

As we have mentioned in the Introduction, one of the 
practically important situations is that of a long charged 
polymer forming complexes with oppositely charged par- 
ticles. Motivated by the problem of nucleosome charge 
inversion, Mateescu et al (1999), Park et al (1999), Sens 
and Gurovich (1999), and Netz and Joanny (1999b) con- 
sidered the complex of a positive sphere with charge q and 
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a negative polyelectrolyte, such as a DNA double helix, 
which has to make some nt > 1 turns around the sphere 
to neutralize it (Fig. 9). These authors predicted a sub- 
stantial charge inversion: more of the polyelectrolyte is 
wound than is necessary to neutralize a sphere. Further- 
more, Mateescu et al (1999) found that a tightly coiled 
polyelectrolyte conformation becomes unstable when the 
chain length exceeds a certain threshold, and then an 
almost straight tail abruptly stretches out (Fig. 9). 

Nguyen and Shklovskii (2001a) emphasized the role 
of correlations in this case of charge inversion. Indeed, 
neighboring turns repel each other and form an almost 
equidistant solenoid, which locally resembles SCL. The 
tail of the polyelectrolyte repels the already adsorbed 
part of the polyelectrolyte and creates a correlation hole, 
which attracts the tail back to the surface (compare Fig. 
6). As a result, the net charge of the sphere with wrapped 
polyelectrolyte, q* , is negative. It is shown that at 
r s — > oo the charge inversion ratio scales as |g*|/g ~ l/n t . 
On the other hand, at small enough r s it can exceed 100%. 




FIG. 9. A polyelectrolyte molecule winding around a spher- 
ical macroion. Due to the Coulomb repulsion, neighboring 
turns, which play the role of Z-ions, are strongly correlated. 

Even more interesting is the system in which charged 
polymer is so long that it forms complexes with many op- 
positely charged particles (see Fig. 4). Examples include 
micelles (Wang et al, 1999), globular proteins (Kabanov 
et al, 1976; Xia and Dubin, 1994), colloids (Braun et al, 
1998; Keren et al, 2001), or dendrimers (Kabanov et al, 
2000; Evans et al, 2001), and, last but not least, his- 
tone octamers forming 10 nm chromatin fiber with DNA 
(Fig. 3). To be specific, we remain with the signs con- 
sistent with the DNA case and consider a long negative 
polymer chain in the solution of positive spheres. If the 
concentration of spheres is large, it is favorable to ad- 
sorb many spheres on the polyelectrolyte chain. As a 
result, each sphere is under-screened by polyelectrolyte 
and has positive net charge. Then, adsorbed spheres re- 
pel each other and the complex forms a periodic necklace 
(see Fig. 4). This necklace is, in fact, a one-dimensional 
Wigner crystal, or SCL, of spheres, which serve as Z- 
ions . Indeed, since the segment of the polyelectrolyte 
wound around one sphere interacts almost exclusively 
with this sphere, it plays the role of the Wigncr-Seitz 
cell. Because of correlations, spheres bind to polyclcc- 
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trolyte in such a large number that the net charge of the 
polyelcctrolytc molecule becomes positive (Nguyen and 
Shklovskii, 2001b). In this case, the charge inversion ra- 
tio scales as Q*/Q ~ n]^ 4 in the absence of monovalent 
salt, where — Q and Q* are the bare and net charges of 
the polyelectrolyte molecule, respectively. This means 
that charge inversion may be larger than 100%. As we 
discussed in Sec. V, charge inversion can be further en- 
hanced by a monovalent salt, in which case Q* /Q ~ n t . 
We shall return to the complexes of a charged chain with 
spheres in Sec. IX. 



VIII. MULTILAYER ADSORPTION 

So far we have not considered the possibility that Z- 
ions fully cover the macroion surface. This may some- 
times happen, particularly when a macroion is very 
strongly charged, or Z-ions are large. Suppose, for in- 
stance, that each Z-ion has a hard core of radius a. In 
this case, the effect of excluded volume of the hard cores 
of the Z-ions adds positive contributions to the surface 
pressure and chemical potential of a SCL (10) that are 
proportional to fc^T and that diverge at the full cover- 
age. Close to the full layer this term compensates and 
then over-compensates the negative Coulomb term /iwCi 
so that charge inversion disappears. Indeed, a full layer 
is incompressible (see Fig. 10b), and, unlike a partially 
filled layer (see Fig. 6 or Fig. 10a), it does not allow for 
the creation of an image-like correlation hole. 

At even larger macroion charge, the second layer starts 
to form, launching a new wave of charge inversion. In 
the beginning, charge inversion is small because all the 
attraction of a new Z-ion approaching the surface is pro- 
vided by a weak interaction with an inflated image in 
the emerging second layer, where once again A^$> a (Fig. 
10c). Continuing, Nguyen and Shklovskii (2001c) arrived 
at the prediction of an oscillating inverted charge Q* as 
a function of Q (see Fig. 10), where charge inversion 
vanishes every time the top layer of Z-ions is full. 

Another way to look at this phenomenon is to examine 
a metallic electrode screened by Z-ions, when the poten- 
tial of the electrode is controlled instead of its charge. In 
this case, oscillations of charge inversion and compress- 
ibility lead to oscillations of capacitance of this electrode 
with the number of adsorbed layers of Z-ions. This is 
similar to oscillations of compressibility and magneto- 
capacitance in the quantum Hall effect, which are re- 
lated to consecutive filling of Landau levels (Efros, 1988; 
Kravchenko et al, 1990; Eisenstein, 1992). In this sense, 
we deal with a classical analog of the quantum Hall effect. 

To conclude this section, let us return to the adsorp- 
tion of weakly charged polyelectrolytes that we discussed 
briefly in the Sec. VII. Dobrynin et al (2001) have shown 
that parallel chains of adsorbed blobs start touching each 



other above the same threshold surface charge density 
a e = ef /I 2 , which corresponds to the onset of squashing 
blobs on the surface. As a result, these authors arrived 
at the conclusion that if weakly-charged chains are ad- 
sorbed on the surface with a > a e , they form a polymer 
liquid. In this liquid, correlations and image formation 
are only due to the uppermost layer, with the thickness 
about that of an unperturbed blob. There are no oscilla- 
tions of inverted charge; instead, charge inversion satu- 
rates at about one layer of blobs and remains unchanged 
afterwards. 




Q/Q 



FIG. 10. Inverted charge Q* as a function of the absolute 
value Q of the bare charge; Qo is the charge of one full layer 
of Z-ions. The dashed line corresponds to the case of Z-ions 
with vanishing radius, Eq. (14). (a) The first layer is not full, 
as in Fig. 6. An approaching new ion creates a correlation 
hole and is attracted to it. (b) The layer is full, there is no 
place for a correlation hole, (c) More than one layer is full. 
A correlation hole exists in the top layer only. 



IX. CORRELATION-INDUCED ATTRACTION OF LIKE 
CHARGES 

The idea of a single screened macroion is a useful one 
in a theoretical context, but in practice it is rarely true 
that there is only one macroion. Typically, there is a cer- 
tain concentration of them, so that interactions between 
them can be important. Let us start with the simplest 
question: consider two macroions, and suppose the con- 
centration of Z-ions in solution is equal to N (see Eq. 
(11)), such that each macroion forms a neutral complex 
with Z-ions. How do these two neutral complexes inter- 
act? It turns out that they attract each other at short dis- 
tances and, therefore, tend to coagulate. In other words, 
two macroions of the same charge may attract each other 
because of the presence of Z-ions. In general, this attrac- 
tion of like charges is as interesting a manifestation of 
correlations as charge inversion, even though our present 
Colloquium emphasizes the phenomenon of charge inver- 
sion. Nevertheless, we must discuss attraction at least 
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briefly, in order to prepare the ground for the subsequent 
discussion of experiments (Sec. X). 

Medium-induced attraction of like charges is nothing 
new in physics, with Cooper pairs of electrons being the 
most prominent example. In the context of molecular 
physics, the most popular explanation of attraction is 
in terms of salt bridges: a divalent ion, such as Mg 2+ , 
can form ionic bonds with two groups with charges of 
— 1 each, in effect connecting them together. This idea 
is indeed adequate if we have, say, two macroion sur- 
faces with regularly placed charges of —1, and there are 
ions of charge +2 between them. However appealing, 
the bridge concept becomes increasingly fuzzy when Z- 
ions have charges of 3 or higher and when charges in the 
macroion are not positioned regularly. 




log M 



FIG. 11. Generic phase diagram of reentrant condensation 
and charge inversion in terms of macroion concentration M 
and Z-lon concentration N (Nguyen and Shklovskii, 2001d). 
Isoelectric composition is shown by the dotted line. The 
dashed "neutrality line" corresponds to neutral complexes in 
the dilute phase. The segregation region is shaded. Minus 
and plus indicate the signs of complexes of DNA with Z-ions. 

Motivated by experimental observations of DNA con- 
densation (see Sec. X), there was a significant effort 
by theorists to try to explain attractive forces by go- 
ing beyond the bridge model. For simplicity, and fol- 
lowing the majority of works, let us consider two planar 
macroion surfaces with some Z-ions between them. Of 
course, Poisson-Boltzmann theory predicts pure repul- 
sion for such a system. However, attraction was observed 
in several computer experiments, including Guldbrand 
ct al (1984), Kjcllandcr and Marcelja (1985), Gronbech- 
Jensen et al (1997), Linse and Lobaskin (1998), and Mor- 
eira and Netz (2000b). On the theory side, several groups 
attempted to go beyond the mean field approximation. 
An important observation is that due to dynamic fluc- 
tuations of counterions, there is an attractive component 
(similar to Van der Waals interactions), but at effectively 
high temperature, or small V (see Eq. for the definition 
of (8)), the Poisson-Boltzmann repulsion still dominates 
and the force remains mainly repulsive (Oosawa, 1968; 
Lau and Pincus, 1998; Ha and Liu, 1998; Podgornik and 
Parsegian, 1998; Golestanian et al, 1999; Golestanian and 



Kardar, 1999). 

On the other hand, attraction of like charges domi- 
nates at effectively low temperatures, when T ^> 1, and 
the idea of spatial correlations between Z-ions, which is 
the central idea of this Colloquium, sheds light on the 
nature of this attraction. Indeed, for extremely large T, 
we deal with two Wigner crystals on the two opposing 
plates; they gain energy when approach each other by 
properly positioning themselves in the lateral direction. 
This was shown by Rouzina and Bloomfield (1996) (see 
also Gronbech-Jensen et al, 1997; Levin et al, 1999; Mor- 
eira and Netz, 2000b). Furthermore, Gronbech- Jensen 
ct al (1997) and Shklovskii (1999a) pointed out that the 
long range order of a Wigner crystal is not important 
for this attractive force. As in the case of charge inver- 
sion, what is important is correlation and short range 
order. As we know, Z-ions form a SCL on the macroion 
surface as soon as Y becomes large. Imagine now bring- 
ing two planar surfaces, along with their respective SCL, 
very close to each other. What we have now is essen- 
tially two copies of Fig. 5 superimposed on top of one 
another, with Z-ions confined in between. Clearly, two 
SCL merge, lowering the energy per Z-ion from e(n) to 
e(2n) < e(n) (see Eq. (9)). Physically, every Z-ion in 
the merged SCL is sandwiched between two macroion 
surfaces, and its Wigner-Seitz cell can be approximated 
as a pair of discs, one on each surface (see again Fig 5). 
The charge of the cell must still be — Ze, but since there 
are two surfaces, the radius of the cell is reduced by the 
factor 1/V2, leading to the energy gain. In some sense, 
this theory returns us to the idea of bridges, albeit on a 
completely new level, with each Z-ion bridging between 
two sides of its Wigner-Seitz cell, which can include many 
surface charges. 

These arguments hold, at least qualitatively, not only 
for plates, but also for macroions of other shapes, in- 
cluding DNA double helices. To be specific, consider two 
DNA double helices. When the concentration of Z-ions 
is equal to N , each DNA is neutralized by Z-ions, and 
the two neutral complexes attract each other at short 
distances. What happens if the concentration of Z-ions 
is higher or lower than A*o? In this case, correlation- 
induced attraction, which is short-ranged, competes with 
Coulomb repulsion, which is much longer-ranged. Note 
that the Coulomb repulsion force is present both at 
N < A^o, when the DNAs are partially screened by Z-ions 
and negative, and at N > A^, when they are overcharged 
and positive. 

What are the implications of this competition between 
attraction and repulsion? They are summarized in Fig. 
11 which shows a phase diagram of the solution with 
number concentrations of macroions, M, and Z-ions, N 
(along with the neutralizing amount of monovalent ions 
and salt.) The major feature of the phase diagram is 
the segregation region, which is the shaded area in the 
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Fig. 11. As the figure indicates, the generic scenario is 
that of reentrant condensation. DNA molecules stay in 
solution and remain negative at N < N C (M), forming 
undercharged complexes with Z-ions. At some concen- 
tration of Z-ions, N — N C (M), repulsion loses to the 
correlation attraction and a condensed phase of DNA is 
formed, coexisting with a dilute phase. The condensed 
phase for DNA represents a (nematic) bundle of helices, 
it exists in the interval N C {M) < N < N d {M). Finally, at 
N = Nd{M), repulsion overcomes the correlation attrac- 
tion, and the DNA molecules dissolve and form positive 
(overcharged) complexes with Z-ions. 

Inside the coexistence region, there is a neutrality line, 
on which the equilibrium dilute phase consists of neutral 
complexes. At very small DNA concentrations, the neu- 
trality condition corresponds to the concentration Ao of 
Z-ions (Eq. 11). To see what happens at larger DNA 
concentration, consider increasing M, starting from the 
overcharged complexes, well above segregation on the 
phase diagram in Fig. 11. When M grows, the solution 
runs out of Z-ions when it approaches the "isoelectric 
line" -TfLM + ZeN = 0, with -77 < and L being the 
DNA linear charge density and length, respectively Near 
this line, the charge of complexes flips sign. Thus, the 
neutrality line crosses over from N = Nq to the isoelec- 
tric line. The border lines N C (M) and Nd(M) follow a 
similar pattern. Although not plotted in Fig. 11, at ex- 
tremely small values of M these two lines join together 
at a critical point, and for smaller M only intramolecu- 
lar condensation of DNA (the coil-globule transition) is 
possible if the DNA molecule is long enough. 

We have considered the phase diagram in Fig. 11 for a 
solution of DNA chains with small Z-ions. As a matter of 
fact, the diagram is qualitatively quite general (Nguyen 
and Shklovskii, 2001d). For instance, it applies to a so- 
lution of DNA with large positively charged particles. In 
Sec. VII, we considered the case of a small DNA con- 
centration, M, and a large concentration of spheres, N, 
which corresponds to the region above the coexistence 
region on the phase diagram in Fig. 11. We found 
that complexes have the form of necklaces, as shown 
in Fig. 4, and that they are overcharged, i.e., contain 
more spheres than necessary to neutralize DNA molecule. 
Large spheres are so strongly bound to DNA that the 
concentration Nq for them is extremely small and any 
real experiment deals with the narrow upper-right part 
of the diagram. Suppose now that there are relatively 
few spheres in the solution, so we are below the neutral- 
ity line. In this situation, chains make an overcharging 
number of turns around each sphere. This is energeti- 
cally favorable due to the repulsive correlations between 
subsequent turns on a sphere surface. The inverted net 
charge of each sphere is about as large as for the case of 
a single sphere discussed in Sec. VII. Furthermore, the 
inverted charge of spheres determines their distribution 



along the chain of polyelectrolyte. Negative spheres repel 
each other and, therefore, the complex once again has 
the periodic beads-on-a-string structure, Fig. 4, which 
resembles the 10 nm chromatin fiber. In the narrow 
vicinity of the neutrality line, even a small correlation 
attraction between touching spheres is sufficient to drive 
aggregation (or coil-globule collapse for long chains) of 
DNA with spheres. 

X. EXPERIMENTAL EVIDENCE OF CHARGE 
INVERSION 

How does the theory of correlated screening compare 
with experiment? For the purposes of this Colloquium, 
we restrict ourselves to a qualitative comparison only. 



4- 




aoi o!i 1 

Polymer Dosage, ^g/mL 

FIG. 12. Mobility of positive latex particles (macroions) in 
the presence of 0.005 M (moles per liter) concentration of 
NaCl as a function of polymer mass concentration. Polymers 
(single-stranded DNA chains) play the role of Z-ions. Differ- 
ent symbols correspond to DNA of the following lengths (in 
monomers): A - 8, o - 10, A - 40, • - 80, □ - 1400. The line 
is drawn to guide the eye. The arrow indicates the isoelec- 
tric point, the polymer mass concentration of 0.17 fig mD 1 
at which DNA charge neutralizes the latex particles. 

First observations of charge inversion were reported 
a long time ago by De Jong (1949). He was able to 
measure electrophoresis of macroscopic aggregates in the 
phase segregation region of the phase diagram Fig 11, 
and observed the reverse of their mobility upon crossing 
the isoelectric line. More recently, Kabanov and his co- 
workers (Kabanov and Kabanov, 1995; Kabanov et al, 
1996; Kabanov and Kabanov, 1998) examined mixtures 
of positive and negative polymers and directly observed 
the interpolyelectrolyte complexes in which a larger poly- 
mer of one sign (playing the role of a macroion) was seen 
to bind an overcharging amount of smaller polymers of 
the opposite sign (playing the role of Z-ions). The ef- 
fect was directly seen due to the reversal of the elec- 
trophoretic mobility of complexes. Furthermore, Wang ct 
al (1999) observed similar reversal for a mixture of poly- 
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electrolytes (macroions) and micelles (Z-ions). Gotting 
ct al (1999) found the reversed mobility for the nanopar- 
ticles (macroions) and short single stranded DNA (elon- 
gated Z-ions). Walker and Grant (1996) demonstrated 
this phenomenon for 120 nm latex particles (macroions) 
with single stranded DNA (Z-ions) ranging from 8 to 
1400 nucleotides; their data are presented in Fig. 12. 
Kabanov et al (2000) and Evans et al (2001) observed 
the reversed electrophoretic mobility for DNA with the 
dcndrimers as Z-ions. 

An interesting observation, apparent from Fig. 12, is 
that the data for different DNA lengths collapse onto a 
single master curve upon re-scaling (in which mobility 
is plotted against mass concentration of DNA instead 
of number concentration, N). We note in passing that 
this observation can be rationalized by the following ar- 
gument. Isoelectric point is obviously determined by the 
number of charged groups per unit volume, but this quan- 
tity is insensitive to the overall length of DNA and is 
simply proportional to the mass concentration. 

All of the above mentioned experimental works rely 
on the reverse of electrophoretic mobility as an indica- 
tion of charge inversion. Indeed, this is conceptually the 
most straightforward approach. It is valid because Z- 
ions are strongly bound to the macroion, with energies 
larger than ksT, and move together with it. On the 
other hand, monovalent ions screening the net charge Q* 
are attracted to macroion with energy much smaller than 
ksT and move in electric field in the opposite direction. 
For these reasons, it is the net charge Q* that determines 
both magnitude and sign of the observed electrophoretic 
mobility. This remains correct also in the case when 
monovalent ions adsorb on Z-ions, effectively reducing 
them to Z'-ions with Z' < Z, as in Fig. 7. In all cases, 
the net charge includes all ions bound with energies in 
excess of fc^T. In a recent molecular dynamics simula- 
tion, Tanaka and Grosberg (2001b) have directly exam- 
ined the mobility of charge inverted macroion complexes 
similar to the one shown in Fig. 7. They confirmed that 
adsorbed Z-ions drift together with the macroion in a 
weak electric field and that the sign of the net charge Q* 
determines the direction of electrophoresis. 

It is also worth saying that the interpretation of elec- 
trophoretic experiments on charge inversion is not af- 
fected by the recent discoveries of Long et al (1996 and 
1998). These authors noted that electrophoretic mobil- 
ity, under some circumstances, may not be entirely deter- 
mined by the charge. For instance, even an overall neu- 
tral object may move in an electric field provided that it 
has strong asymmetry of charge distribution. A simplest 
example consists of two balls of different radiuses, rigidly 
connected by a thin rod, and having opposite charges. 
The effect is due to the fact that external electric field 
acts not only on macroion itself, but also on surround- 
ing co- and counterions, causing the latter to flow and 



to exert viscous friction drug forces on the macroion. 
Strong geometrical asymmetry of positive and negative 
charges leaves these drug forces unbalanced. In this col- 
loquium, we discuss only macroions in which the bare 
charge is uniformly distributed on the surface. Then, the 
charge of Z-ions is practically uniform, too (apart from 
the small length scales of the order of R, the distance be- 
tween neighboring adsorbed Z-ions). For such macroions, 
reversal of electrophoretic mobility does indeed indicate 
the inversion of charge. 

For a more detailed comparison with experiments, we 
should remember that charge reversal is expected to be 
accompanied by coagulation, as discussed above in Sec. 
IX. Whether in equilibrium or not quite in equilibrium, 
these large complexes should scatter light strongly. There 
are many experiments reporting such observations. 

Let us begin with DNA. It has been known for some 
time that at some critical concentration of Z-ions, N c , 
DNA abruptly condenses into large bundles (Bloom- 
field, 1996). Recently it was discovered that at a much 
larger critical concentration, N d , the bundles dissolve 
back (Saminathan et al, 1999; Pelta et al, 1996a; Pelta 
et al, 1996b; Raspaud et al, 1998; Raspaud et al, 1999). 
Specifically, for the spermine ions (Z = 4), it was found 
experimentally that N c = 0.025 mM and N d = 150 mM. 
If one interprets these data in the framework of the the- 
ory of Nguyen et al (2000c), the experimental values of 
N c and Nd imply that, for spermine, N — 3.2 mM and 
the binding energy of two helices per one spermine ion is 
u = Q.'iksT. The last value agrees with the one obtained 
by a different method (Rau and Parsegian, 1992). 

Let us now discuss some other systems. Wang et al 
(1999) studied complex formation in mixture of micelles 
and oppositely-charged polyelectrolytes. In this exper- 
iment, the total charge of micelles was controlled by 
changing the concentration of the cationic lipid in the 
solution. In agreement with the above theory, measure- 
ments of dynamic light scattering and turbidity (coeffi- 
cient of light scattering) show that complexes condense 
in bundles and the solution coacervates in the vicinity of 
the point where mobility crosses over between two almost 
constant values, positive and negative. 

For the complexes of latex particles with DNA of vari- 
ous lengths, examined by Walker and Grant (1996), equi- 
librium conditions were not found, but a significant rate 
of aggregation of latex particles was observed in the same 
narrow range of DNA concentrations where the mobility 
flips the sign. 

There is a large body of interesting experimental 
(Radler et al, 1997; Koltover et al, 1999) and theoret- 
ical (Harries et al, 1998; Bruinsma, 1998) work on phase 
diagrams and overcharging of lamellar cationic lipid- 
DNA self-assembled complexes. These solutions always 
seem to have aggregates due to hydrophobic attraction 
of lipids. A phase diagram of this kind has been sketched 
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by Radler (2000). 

Another large group of works with nice (and tech- 
nologically useful) examples of overcharging deals with 
the sequential adsorption of multiple layers of polyelec- 
trolytes of alternating sign (see Decher (1997) and ref- 
erences therein). Say, a positive surface is first treated 
with negative polymers, gets overcharged and becomes 
negative, then it is treated with positive polymers, gets 
overcharged etc. This procedure works reliably with up 
to many tens of layers. Theoretical interpretation of 
this technique was discussed by Castelnovo and Joanny 
(2000). It falls beyond the framework of the present Col- 
loquium as it involves certain kinetic considerations while 
we deal in this Colloquium only with equilibrium phe- 
nomena. 

XI. CORRELATIONS "IN A SHEEP'S SKIN" 

We mentioned correlations so many times in this Col- 
loquium, that the reader may want to ask: Are there 
alternative, correlation-independent, electrostatic mech- 
anisms leading to charge inversion? Our answer is no, 
and we argue that the correlation-based mechanism is 
the universal one. This fact notwithstanding, we should 
say that correlations may show up in a number of ways, 
sometimes masked like a wolf in a sheep's skin. 

To understand better the role of correlations, let us 
first consider the case of no correlations. Namely, sup- 
pose, we have a set of randomly positioned point-like 
charges, equal numbers of +e and — e. It is easy to es- 
tablish that the averaged interaction energy in such sys- 
tem is exactly zero, where average is taken over indepen- 
dent random positions of all charges. Similarly, averaged 
electric field in the system is also zero. Both these state- 
ments remain correct also for the one component plasma, 
in which point-like charges of one sign are randomly po- 
sitioned on the smeared uniform background of the op- 
posite sign. 

The essence of screening is that charges in plasma are 
not positioned randomly. Correlations happen because 
ions reconfigure themselves non-randomly to gain some 
energy. This so-called correlation energy is well known 
in plasma physics (Landau and Lifshitz, 1977). It is neg- 
ative, meaning that correlated configuration is lower in 
energy, or more thermodynamically favorable, than ran- 
dom configuration. 

These simple facts allow us to understand the under- 
lying role of correlations in one of the theories suggested 
in the literature to explain charge inversion. This the- 
ory, put forward by Park et al (1999), views monovalent 
counterion release from Z-ions as the driving force behind 
charge inversion. We argue that while counterion release 
accompanies, it is itself driven by correlations. As we 
did already more than once, let us imagine that DNA 



molecules, along with their Onsager-Manning-condcnsed 
small ions, are being adsorbed on the macroion one at a 
time (possibly releasing some of their countcrions). Let 
us further consider the moment when the neutralization 
condition has just been achieved. Pretend now that DNA 
rods (Z-ions) are distributed randomly on the surface, 
uncorrelated in both positions and orientations. In this 
case, next arriving DNA molecule feels no average field, 
so that it has no reason to release its counterions. The sit- 
uation is completely different if DNA molecules are corre- 
lated on the surface (see Figs. 2 or 8), where locally each 
molecule is surrounded by a correlation hole - the pos- 
itive stripe of the background charge (the Wigner-Seitz 
cell). The corresponding field, or the positive potential 
of the Wigner-Seitz cell, causes the release of counterions 
from DNA not only at the neutrality point, but even if 
the surface overall is overcharged (see the solution of this 
problem given by Eqs. (20), (21)). These qualitative 
arguments can be formulated also more quantitatively 
(Nguyen et al, 2001). Thus, correlation hole, or adjust- 
ment of DNA molecules to each other, or image charge, or 
correlations (all synonyms!) is the driving force for both 
counterion release and charge inversion; under a sheep's 
skin of countcrions release, there is a wolf's face of cor- 
relations. 

Let us now discuss another approach which we call 
metallization. It was pioneered by Mateescu et al (1999), 
who considered complexation of a polyelectrolyte with 
a sphere, and by Joanny (1999), who examined adsorp- 
tion of flexible polymers on a charged plane. Metalliza- 
tion theory considers adsorbed Z-ions as a continuous 
medium similar to a metal, while still treating the bulk 
solution as consisting of discrete charges. 

We argue that metallization theory is in fact an over- 
estimate of correlation effects. Indeed, as we illustrated 
with Fig. 6, SCL behaves as a metal on the length scale 
above R, the distance between neighboring adsorbed Z- 
ions. It behaves as a metal in the sense that it responds 
to the approaching new ion by forming an image. Clearly, 
the smeared continuum is also a metal in this sense, and 
even a better metal - it is a metal on all length scales. 
Another way to view that same physics is to note that 
correlation suppresses electric field of every Z-ion beyond 
certain distance of order R, while for the smeared con- 
tinuum the field is suppressed everywhere and, therefore, 
the effect is overestimated. 

The latter view suggests also another fruitful interpre- 
tation of metallization approach, which is in terms of the 
sclf-cncrgy of Z-ions. Here, we resort to the terminol- 
ogy in which self-energy of an ion is identified with the 
energy of the electric field outside certain cutoff length, 
such as an ion size. (We are not interested in the electric 
field on the smaller length scales: although energy of this 
field is infinite, it does not change upon any processes 
considered here, such as adsorption of Z-ions.) Using 
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this language, the metallization theory becomes physi- 
cally transparent. Indeed, in this theory adsorption of 
Z-ions is energetically favorable because they do have 
self-energy while in the bulk solution and loose it com- 
pletely upon adsorption. Once again, real screening, or 
correlations between Z-ions, correspond to suppression 
of a part of self-energy, corresponding to the energy of 
the electric field beyond the distance of order R. Energy 
of the field between the ion size and R is the amount of 
overestimation by the metallization theory. By contrast, 
the Poisson-Boltzmann approximation fails to describe 
charge inversion precisely because it smears Z-ions ev- 
erywhere and neglects their self-energies. 

Representation of correlations in terms of self-energy 
allows us to address one more practically important prob- 
lem, namely, the discreteness of charges (Nguyen and 
Shklovskii, 2001e). Indeed, instead of uniformly charged 
surface, like Fig. 1, it would be closer to reality (and 
to a chemist's heart) to think of a macroion as having 
some charged groups with unitary charge each. For in- 
stance, negative (— e) charges of DNA are positioned on 
an external spiral rim of the double helix, at the distance 
A = 0.67 nm from each other along the rim. 




FIG. 13. Schematic representation of charge fractionaliza- 
tion. a) One strand of negative charges of DNA (empty balls) 
is completely neutralized by positive Z-ions with Z = 3, their 
charged groups are shown by black balls. A new Z-ion is ap- 
proaching, b) The new Z-ion is "digested." Its charge is split 
in +e charges of Z defects, tails and archs. 

To be specific, let us consider a macroion which is a 
regular lattice of charges — e ("unfolded DNA strand") 
and Z-ion having also a linear array of charges +e (a 
short polyelectrolyte). Importantly, Z-ions always have 
some degree of flexibility; for instance, in the case Z = 3, 
as shown in Fig. 13, we can imagine that they freely 
bend in the middle. For simplicity we assume that the 
distance between charges in the Z-ion, B, matches ex- 
actly that in the macroion, A: A = B. Finally, we as- 
sume that the Z-ion charges and the macroion charges 
can approach each other to the minimal distance much 
smaller than A. Then Z-ions can attach to the macroion 
locally compensating each charge and, therefore, achiev- 
ing complete neutralization, as we show in Fig. 13a. The 
neutralization is so perfect that it is difficult to imagine 



how another Z-ion can be attached. Fig. 13 explains 
why this happens. Similarly to Fig. 6, suppose that the 
macroion is already neutralized by Z-ions and a new Z- 
ion comes. Then it turns out energetically favorable to 
disturb order among the neutralizing Z-ions by Z defects 
in Z independent places, thus opening a room for a new 
Z-ion. In each defect, one charge of Z-ion is detached 
from the corresponding macroion charge, forming posi- 
tive tail or arch above the surface and leaving negative 
vacancy on the macroion. Then, shifting Z-ions along the 
macroion, Z vacancies can join together and form a large 
vacancy capable to accommodate an entire new Z-ion. A 
net result is that Z disconnected charges +e appear on 
top of completely neutralized macroion (Fig. 13b), or, 
in other words, the charge of Z-ion is fractionalized. To 
avoid misunderstanding, we emphasize that none of the 
chemical bonds is really cut. 

Fractionalization effectively eliminates self-energy of 
the free Z-ion. Indeed, self-energy of the Z-ion is simply 
the energy of repulsion between Z positive constituent 
charged groups in the extended conformation which the 
Z-ion assumes in the solution. In the fractionalized state, 
charges are far apart and practically do not interact, so 
that self-energy is gained. These results are easily gener- 
alized for the case B < A, when Z-ions have larger linear 
charge density than macroion. For instance, if B = A/2, 
then Z-ions repelling each other form a " Wigner crystal" 
on top of the lattice of macroion charges where they al- 
ternate with vacant places (similarly to Figs. 1, 5, 6, and 
8). 

Is fractionalization a correlation-independent mecha- 
nism of charge inversion? Of course, not: this phe- 
nomenon is solely due to the correlated distribution of 
Z-ions, which avoid each other at the macroion. Frac- 
tionalization is yet another mask under which correla- 
tions may show up. 

XII. CHARGE INVERSION IN A BROADER PHYSICS 
CONTEXT 

In conclusion, we would like to show that charge in- 
version studied in this Colloquium has many physical 
analogies. There are other "charge inverted" systems 
in physics. Let us start from the hydrogen atom. It 
is known that it can bind a second electron, forming the 
negative ion H _ with an ionization energy of approxi- 
mately 0.05Ry (Massey, 1938) . We can consider this 
effect as the inversion of proton charge. Attraction of 
the second electron to the neutral atom is due to the 
Coulomb correlation between electrons: the first electron 
avoids the second one, spending more time on the oppo- 
site side of the proton. In other terms one can say that 
binding is related to polarization of the neutral core. 

Negative ions - nuclei overcharged by electrons - ex- 
ist also for larger atoms. Mean field Thomas- Fermi or 
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Hartree theories fail to explain negative ions (Landau and 
Lifshitz, 1980). One must include exchange and Coulomb 
correlation holes to arrive at a satisfactory theory ex- 
plaining bound state and the nonzero ionization energy 
of a negative ion (Massey, 1938). The Thomas- Fermi 
theory of an atom is an analog of the Poisson-Boltzmann 
theory of electrolytes. It is not surprising that both fail 
to explain charge inverted states. 

Similar considerations apply also for a macroscopic 
metallic particle. Electrons in such a particle have a 
negative (compare to vacuum) chemical potential or, in 
other words, positive work function. The work function 
is known to vanish in the Thomas-Fermi or Hartree ap- 
proximations (Lang, 1973). 

The energy |/xwc| plays the same role for Z-ions on 
the insulating macroion surface as the ionization energy 
of a negative ion or the work function of metallic parti- 
cle for electrons. Similarly to electrons, charge inversion 
of a charged insulating macroion by Z-ions can not be 
obtained in the mean field Poisson-Boltzmann approxi- 
mation. Only correlations of Z-ions on the surface of the 
macroion can lead to charge inversion. 

Let us now return to the Onsager-Manning conden- 
sation (Manning, 1969; Sec. II). Kosterlitz and Thou- 
less (1972) discovered a similar threshold phenomenon 
for generation of free vortexes in two dimensional super- 
fluids or superconductors. They noticed that due to the 
logarithmic form of attractive interaction, two vortexes 
of opposite sign decouple only above some critical tem- 
perature, Tkt- Later, Kosterlitz- Thouless theory was 
applied to unbinding of dislocations and disclinations in 
the theory of defect-induced melting of two-dimensional 
crystals (Nelson and Halperin, 1979; Young, 1979). 

In the Kosterlitz- Thouless theory one can identify ana- 
log of the short-range correlation contribution to the 
chemical potential of Z-ions (which we called |/Ltvpc|)- 
This is the energy of creation of the two vortex cores. 
Similarly to |^nyc| 7 this energy provides additional bind- 
ing of vortexes and strongly reduces the concentration of 
free vortexes at T > Tkt (Minnhagen, 1987). In con- 
trast to the Kosterlitz-Thoulcss theory, the short range 
contribution |/iyi/cl was only recently introduced (Pcrcl 
and Shklovskii, 1999). 

Another physical analogy already mentioned in Sec. 
VIII is the integer quantum Hall effect. We should add 
that charge fractionalization, as illustrated by Fig. 13, is 
an analog of the fractional quantum Hall effect (Prange 
and Girvin, 1990). Finally, Fig. 13 bears also analogy 
with electron charge fractionalization in the polyacetilcne 
(Brazovskii and Kirova, 1991). 

XIII. CONCLUSIONS AND OUTLOOK 

To conclude, we have discussed the physical picture 
of screening for the case of strongly interacting ions. 



This case appears to have been overlooked for many 
decades, since Debye. Its theoretical study, motivated 
by the experiments mainly on gene delivery and chro- 
matin structure, revealed new interesting physical in- 
sights. Specifically, correlations between screening ions 
lead to such counterintuitive phenomena as charge inver- 
sion, reverse electrophoretic mobility, attraction of like 
charged molecules or colloids, etc. The physical theory of 
these phenomena is aesthetically attractive as it presents 
many parallels with other areas of physics, ranging from 
quantum Hall effect to atomic physics and metals. The 
potential applications are many, in both chemical and 
biological realms. We can mention here all sorts of ma- 
nipulations with DNA, both for biological purposes and 
for using DNA as an assembly tool for non-biological nan- 
otechnology. They all require understanding of electro- 
static properties of DNA chains. Many diseases have to 
do with mis-assemblies of charged proteins, such as actin; 
we need to understand better the assembly of such ob- 
jects. Food, cosmetic, paper, and waste water treatment 
industries, are all about charged colloids, and the list of 
applications is easy to continue. In brief, this theory is 
one of the busy junctions where physics meets chemistry 
and biology. 
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